Frustrated local moment models for Fe-pnictide magnetism 
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The low energy spin excitations of the Fe pnictide parent compounds have been determined 
by inelastic neutron scattering and interpreted within the local moment Ji a ,b-J2 Heisenberg model 
with orthorhombic symmetry. This has led to alternative exchange models that strongly differ in the 
size of anisotropy. Although the compounds are itinerant the localized spin model can explain basic 
features of the excitations. The inherent frustration of this model leads to quantum fluctuations and 
possible moment reduction. We investigate this question in detail using spin wave approximation and 
partly exact diagonalization Lanczos calculations for finite clusters. We find that the orthorhombic 
anisotropy stabilizes the columnar AF phase and its moment. For the exchange models proposed 
from inelastic neutron scattering we can exclude a strong influence of frustration on the moment 
size. We also investigate dependence of magnetization and susceptibility on field and temperature. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.30.Ds 



I. INTRODUCTION 

The discovery of Fe pnictide superconductors has given 
new impetus to study the interrelation between mag- 
netic and superconducting instabilities in condensed mat- 
ter. Previous investigation of strongly correlated heavy 
fermion and cuprate compounds have also shown this 
connection. However, unlike those compounds the pnic- 
tides are only moderately correlated 1 - and in contrast to 
cuprates already the parent compounds are metallic. At 
lower temperatures they exhibit a structural phase tran- 
sition from tetragonal to orthorhombic symmetry and 
simultaneously or subsequently show magnetic order 2,3 . 
For 122 and 1111 compounds the latter is found to be of 
the columnar AF type in the FeAs planes corresponding 
to wave vector Q = (irfi) which is equal to the nest- 
ing vector of hole and electron Fermi surface pockets ob- 
tained from density functional calculations of the elec- 
tronic structure. This commensurate magnetic structure 
is stable even for a considerable range of dopings 4,5 . How- 
ever it was noticed 4 * 5 - that the size of the staggered mo- 
ment depends strongly on the details of the calculation, 
especially on the out-of-plane position of As atoms, and 
is much larger as the experimental values obtained from 
neutron diffraction which is less than 1/ib per Fe. The 
moments are oriented parallel to the ordering wave vector 
which is aligned to the long (a) axis. 

On the other hand the results of inelastic neutron scat- 
tering (INS) on the the low energy spin excitations have 
shown that they may be successfully described within 
a localized Heisenberg model that includes interactions 
up to next nearest neighbors, i.e., exchange bonds along 
the sides and diagonals of the FeAs layers^—. There- 
fore the localized Heisenberg model is of the anisotropic 
Jia,b-^2 type with different exchange constants along the 
orthorhombic a,b axes. The local moment picture has 
a further interesting aspect: Depending on the ratios 
J\a,bl Ji the magnetism may exhibit frustration which 
can strongly reduce the size of the ordered moment. As 
mentioned above the observed moments are much smaller 



than the calculated ones. Experimentally they vary from 
0.36/xb in LaFeAsO to 1/ib in SrFe2As2 whereas LDA 
calculations give and 1.7/ie 5 respectively if the 

experimental values for the As positions are used. This 
translates to a relative moment reduction factor of 0.19 
to 0.59. Within a model that includes only spin degrees 
of freedom one may conjecture that magnetic frustration 
and associated enhanced quantum fluctuations are the 
source of the strongly reduced ordered moments in the 
Fe pnictides. Such an idea, however needs to be treated 
with care. Firstly frustration itself is not a well defined 
concept for itinerant spins, but rather refers to the local 
moment model. If the latter is indeed used the extent 
of frustration depends crucially on the ratios of exchange 
constants and their anisotropies. 

No consensus on the proper exchange model in the 
local moment picture of Fe pnictides has emerged. Ba- 
sically two proposals, both from INS- - — and local den- 
sity calculations 4,5 have been made: In the first choice 
Jib ~ Ji a , in the second one \Ju\ *C Ji - In principle 
this will lead to very different dispersion for wave vectors 
along the 6*-axis in the two cases which may be checked 
experimentally. But the case remains undecided sofar. 
The second very anisotropic case (Jib was even reported 
slightly ferromagnetic 5 ) is hard to understand in terms 
of a small tetragonal-to orthorhombic structural distor- 
tion with (a — b)/a ~ 0.5 • 10~ 2 . This indicates that 
the exchange constant should not be interpreted as bond 
exchange energies. Rather they are fit parameters ob- 
tained from mapping the total LDA energy with spiral 
magnetic structure to the classical ground state energy 
of a localized Heisenberg model 4 * 5 -. 

One important motivation why a localized model is 
nevertheless worthwhile to study was provided by the 
INS results. In a metal one would naively expect that 
spin waves with larger wave vectors should quickly be- 
come overdamped when they merge with the continuum 
of particle hole excitations. In fact it was found that well 
defined spin waves exist throughout the Brillouin zone 1 ^ 
giving support for the local moment picture. Even in 
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weakly correlated 3d-compounds local moments may be 
stabilized due to Hund's rule coupling in a multi-orbital 
case like Fe pnictides as has been shown in Ref. HH This 
is not unlike the situation in elemental Fe ferromagnetism 
where it was found that Hund's rule exchange stabilizes 
the local momen ta 13 ' 14 . In fact spin waves in elemental 
fee Fe also exist throughout the Brillouin zone and may 
be described by a localized Heisenberg model although 
Fe is a good metali^. The itinerant nature of magnetism 
in the Fe pnictides does not by itself exclude that a local 
moment model is a good starting point for studying the 
low lying spin excitations. Further support for this con- 
jecture was given by functional renormalization group 
(FRG) calculations for the isotropic {J\ a = Ju = J\) 
model 16 . Starting from a multiorbital extended Hubbard 
model it was shown that the J1-J2 model is a valid de- 
scription for the dominant low energy correlations for a 
wide range of parameters. Coexistence models for both 
localized and itinerant moments in the pnictides have 
been proposed in Refs. H3 andfTH 

Extended models including orbital degrees of freedom 
within a localized Kugel-Khomskii type approach have 
been proposed in Refs. [I9I - I22I In this case the ground 
state can exhibit orbital order which may lead to effective 
orthorhombic exchange anisotropy and low ordered mo- 
ment. Furthermore itinerant multi-orbital models 2 ^—, 
also including the effect of orbital orde r 12 ' 26 ' 27 for the 
magnetic ground state have been proposed. 

The INS experiments^— were interpreted within mod- 
els including only spin degrees of freedom. Since we want 
to refer the exchange parameters found there we also re- 
strict to these type of models. In this work we have made 
a systematic survey of the anisotropic (orthorhombic) 
two dimensional J\ a ,b-J2 Heisenberg model. Our main 
subject is to examine carefully to which extent the 'frus- 
tration' of nearest and next nearest exchange constants in 
this model plays a role in the reduction of the staggered 
moment as observed in neutron diffraction. For that pur- 
pose we are using both analytical spin wave calculations 
and numerical Lanczos method for finite clusters of the 
two-dimensional (2D) rectangular lattice. In the central 
part we investigate the evolution of the staggered mo- 
ment, reduced by quantum fluctuations as function of the 
frustration and anisotropy ratios. This allows us to make 
a quantitative evaluation of the importance of frustration 
in the local moment model for Fe pnictides by comparing 
with the results for the experimental exchange constants. 
The isotropic model with J\ a ,b = J\ has been studied pre- 
viously within various approximation s 28 ' 29 and has also 
been extended including interlayer couplin g 30 ' 31 . Spin 
waves for the anisotropic model in zero field are also dis- 
cussed in I32I 

In Sect. |H]we introduce the localized spin model for 
the 2D orthorhombic Fe pnictide layers and discuss its 
parametrization. In Sect. IIIII we calculate the ground 
state energy and phase diagram of the model and the 
corresponding location of known Fe pnictide compounds. 
The reduction of the ordered moment by quantum fluctu- 



ations using spin wave expansion is discussed in detail in 
Sect. lIVl The effect of an external field in the anisotropic 
model is addressed in Sect. El Finally Sect. IVII gives the 
summary and conclusion. 



II. LOCALIZED SPIN MODELS FOR 2D 
RECTANGULAR LATTICE 

We start from a localized spin model with effective spin 
size S = 1/2. The latter is suggested by INS and LDA 
results given in Table U which are compatible with this 
value. A stronger argument is given by the Gutz wilier ap- 
proach to the multiorbital Hubbard model 12 which sug- 
gests that 2S ~ 1 for reasonable model parameters. 

The orthorhombic symmetry allows for different n.n. 
exchange parameters J\ a ,b, however the extreme differ- 
ence for some parameter sets in Table|I]can hardly be jus- 
tified by the simple effect of exchangestriction on nearest- 
neighbor (n.n.) exchange caused by the orthorhombic 
distortion. As mentioned above, a more likely source 
of large exchange anisotropy is the presence of underly- 
ing orbital order which may appear simultaneously with 
magnetic order. However a quantitative prediction of the 
amount of anisotropy seems difficult. Due to the large 
differences in the proposed exchange parameters we treat 
the anisotropy as a free parameter. 

INS results also show the existence of a small spin gap 
of the order ~ lOmeV. For a S = 1/2 system this can 
most easily be modeled by a uniaxial out-of-plane spin- 
space exchange anisotropy which will be included in the 
model for completeness but not discussed in detail. Fur- 
ther insight in the underlying frustrated exchange model 
may be gained from field field dependence of magnetiza- 
tion and susceptibility as has been shown for a different 
class of compounds— Therefore we also include a Zee- 
man term in the model. 

The effective localized spin Hamiltonian we shall dis- 
cuss in this paper then has the form 

(<y> 1 

where the sum in the first term extends over bonds (ij) 
connecting sites i and j. We assume an interaction in 
spin space of the form 

J,, diag (./,;.. /,;../;,). (2) 

To conserve U(l) symmetry, the magnetic field points 
into the z direction defined by the anisotropy introduced 
above. Suppressing the direction index we set 

!Ji a if Rj = Ri±e x 
Ju if Rj = Ri±£ y , (3) 
J 2 if Rj = Ri ±e x ±e y 

i.e., we restrict Eqs. (JTJ and © to nearest- and next- 
nearest neighbor exchange on a rectangular lattice. For 
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the discussion of the complete phase diagram of this 
Hamiltonian, we use a more convenient parameterization 
of the exchange terms and write 

Ji a — V2J C cos (ft cos 9, 
Jib = v2 J c cos (ft sin 9, (4) 
J 2 = J c sin (ft, 

introducing an overall energy scale J c (this should not 
be confused with J z ), a frustration angle (ft, and an 
anisotropy parameter 9. For 8 = ir/A, the above Hamil- 
tonian reduces to the square-lattice case {Jia — Jib) in~ 



vestigated before, e.g. in Refs. [331 - 4351 for V oxides. 



III. CLASSICAL PHASES AND GROUND 
STATE ENERGIES 

For the isotropic model (9 = 7r/4) these are well known 
and serve as starting point for discussing frustration ef- 
fects. We first find out how their energies and existence 
regions are modified for the general anisotropic case with 
— 1 < 9/n < 1. On each site i, we introduce a local co- 
ordinate system, where the z' axis is oriented parallel to 
the local magnetic moment, and we have 



Sf 
Sf 
Sf 



cos(QR t ) -sm(QRi) 0\ /cos© -sin© 
sin(Qi? I ) cos(Qi^) 1 

l / \ sin© cos© 



Sf 
Sf 
Sf 



(5) 



with the ordering vector Q in the xy plane perpendicular 
to the magnetic field which points along the z axis. At 
finite values of H , the spins form an umbrella-like struc- 
ture around the direction of H; the respective canting 
angle (not to be confused with the anisotropy param- 
eter 9) is measured relative to the field direction (global 
z axis); = corresponds to the fully polarized state, 
and = 7r/2 to the state with vanishing magnetic field. 

With h = guBH, the classical Hamiltonian then reads 



Hd = NS 2 



Ji (Q) + A(0) cos 2 - - cos 



(6) 



where we have introduced the Fourier transform 



Uk) = - £ j9, e -m-^) = 1 j2 J>- iUn (7) 



for a = {-L,z} and the last sum runs over all bonds n 
connecting a fixed site i with its neighbors. The coeffi- 



cient A(0) = A(k = 0) is defined by 



A(k) = J z (k) + -[J±(k + Q) + J ± (k-Q) 



2J±(Q). 
(8) 

The reason for defining A(k) in this way will become clear 
later in Sec. IIVI Minimizing Eq. ^ with respect to 0, 
we get the classical canting angle C via 



cos C 



2SA(0) ' 



(9) 



and Eq. © reads 



Ha = NS 2 [j_l(Q) - A(0) cos 2 C . (10) 
Minimizing this Hamiltonian with the exchange param- 
eters from Eq. © with respect to the components of Q 
leads to the four well-known classical phases with order- 
ing vectors 



Q 



ferromagnet (FM) 

(w/a,Tr/b) Neel antiferromagnet (NAF) 

(ir/a, 0) columnar AF along a (CAFa) 

(0, n/b) columnar AF along b (CAFb) 



The minimization condition reduces to dJ±(Q)/dQ = 
and is thus field independent and depends on the trans- 
verse exchange parameters only. The classical ground- 
state energies are 



E± = NS 2 



I 1 - + T ± 

J la + J lb 

2J 2 — (Ji a + Jj^l 

Jib ~ (Jla + ^ J 2 ) 
Jla ~ [Jib + 2J 2 ) 



2 ( J 2 Jf t 

^lb ~ Jib) 
J la J la) 



{Jla + Jla 
[Jla + Jla 
(Jib + Jib 



~ (Jl\ 

2 (Jit 
2UJ- 



Jib) 

-Ji) 
-Ji) 



cos 2 C 

COS 2 C 
COS 2 C 



FM 
NAF 
CAFa 
CAFb 



(12) 
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where cosO c = h/h s , and the critical or saturation fields for the nonuniform phases are given by Eq. 




NAF 

CAFa . (13) 
CAFb 



r 



The minimization condition contains an additional ex- 
tremal solution having an incommensurate wave vector 
given by 



Ju 



, cos Qyb 



2T2'' 



with a ground state energy 



Jla-J 



lb 



8J2 



(14) 



(15) 



However, this energy for the incommensurate wave vector 
is always higher than or equal to the energy in Eq. (|12[) 
of the commensurate ground state corresponding to the 
values chosen for the exchange constants. 

From the classical ground state energy in Eq. (fT2|) 
one may already construct the phase diagram in the <j), 9 
plane, however for the following discussions it is impor- 
tant to include the effect of quantum fluctuations. 



IV. QUANTUM FLUCTUATIONS AND 
ORDERED MOMENT SIZE IN SPIN WAVE 
APPROXIMATION 

In the regions of the <f>, 9 plane where two or more of 
the classical phases become degenerate large quantum 

I 



fluctuations appear and reduce or suppress the ordered 
moment. These are the strongly frustrated regions of 
the phase diagram. One may approach them to some 
extent by starting from the stable region, and calculate 
the contribution of zero point fluctuations in spin wave 
approximation. This leads to an improved ground state 
energy, an estimate for the reduction of the ordered mo- 
ment and for the extent of the instability region where 
magnetic order breaks down. This program has been suc- 
cessfully implemented before for the isotropic cas e 33 ' 34 . 
It will now be carried out for the more general model 
in order to quantify the importance of frustration and 
quantum fluctuations for the compounds listed in Table 

in 



Returning to Eq. (JXJ) expressed in the local coordi- 
nate system introduced in Sec. IIII1 Eq. ([5]), we apply 
a Holstcin-Primakoff transformation and carry out a 1/S 
expansion, keeping terms up to first order in 1/S. (We 
regard h formally as proportional to S.) Next, we ap- 
ply a Fourier transformation. A detailed description of 
the necessary steps can be found in Appendix [A] The 
resulting Hamiltonian is given by 



H=H c i + NS 
S 



■E 



J± (Q) + A(o) cos 2 e - — cos e 



+ V2NS 3 sin8 



A(k) - cos 2 6 \B(k) + 2,4(0) J + - cos 9 



A(0) cos 6 



25 



a l a k + a -j? flt _fe 



+B(k) (1 - cos 2 6) (a % a_ % + a|a f _-) + C(k) cos 6 (ata £ - a_ga^) } 



(a + ajj 



(16) 



where A(k) is defined in Eq. ((SJ), and 

B(k) = J z (k) - 1 (j^k + Q) + J±(k-$j) , (17) 

C(k) = J ± (k + Q)- Jx(fc-Q). (18) 

Eq. (ITB1) still contains a part which is linear in the bosons 
with zero momentum. It occurs only in finite magnetic 
fields. In equilibrium, when 6 = C (see Eq. ([9])), this 
part vanishes, and we will drop it in the following. A 
subsequent Bogoliubov transformation leads to the final 



r 



form 



H 



rtcl + "H zp - 
where H c i is given by Eq. 



S^2E(h,k)al 



TO]), and 

S 
2 



n zp = NSJ ± (Q) + -J2E(h,k) 



(19) 



(20) 



is the zero-point energy contribution to the total ground 
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state energy E gs — Hci+'Hzp- The former corresponds to a ' -. = v^a^ + u^a' -,. (22) 

the magnon excitations described by the boson operators 

The k sums in the equations above span the full crystal- 
lographic Brillouin zone. For the spin-wave dispersion, 
a k = u k a k V k a ^ p (21) we obtain the expression 
I 

2 + C(fc)cos6 c . (23) 



E{h,k) 



A{k) - B(k) cos 2 9 C 



B(k) (1 -cos 2 6 C ) 




C(k) only occurs at finite magnetic fields, and because 
it is antisymmetric in k, it does not contribute to the 
zero-point fluctuations. 



A. Total ground-state energy 

We now calculate the total ground state energy in spin 
wave approximation to determine the zero-field phase di- 
agram. We also will give a comparison to the classical 
ground state energy and the results for finite clusters ob- 
tained from the exact diagonalization Lanczos method. 

Unless explicitly mentioned otherwise, we assume = 
C from here on. Furthermore spin space anisotropy is 
ignored (J^j = J z ?) from now on. The ground-state en- 
ergy is given by the sum of the classical energy, Eq. (fTXJj) 
and the zero-point fluctuations of the magnons, Eq. (|2T))) . 
with the dispersion from Eq. (|23| . Explicitly we have for 
isotropic exchange parameters 

FM: no zero-point fluctuations, 
J{Q) = Jia + Jib + 2J 2 , 

A(k) = 2 [J Xa {cos{k x a) - 1 ) + J u (cos(fc y 6) - 1) 

+2 J 2 (cos(k x a) cos{k v b) - 1)] , (24) 
B(k) = 0, 

NAF: Q = (tt.tt), 

J(Q) = -(J la + J lb ) + 2J 2l 

A(k) = 2 [J la + J lb + 2 J 2 (cos(fc^a) cos(k y b) - 1)} , 
B{k) = -2 [Ji a cos(fc^a) -I- Ju cos(k v b)] , (25) 

CAFa: Q = (tt,0), 

J(Q) = Jib - (Jla + 2J 2 ) , 

A(k) - 2 [J la + J lb (cos(k y b) - 1) + 2 J 2 ] , (26) 
B(k) = —2 [Ji a cos(k x a) + 2 J 2 cos(k x a) cos(k y b)] . 

CAFb: Q = (0,vr), 

J(Q) - Jla ~ (Jlb + 2 J 2 ) , 

A(k) - 2 [J la (cos(k x a) - 1) + J lb + 2J 2 ] , (27) 




FIG. 1. Ground-state energy in linear spinwave approxima- 
tion of the frustrated Heisenberg Hamiltonian on the rectan- 
gular lattice as a function of the frustration angle <f> and the 
anisotropy parameter 9. Energy unit is the overall energy 
scale J c , the magnetic field is zero. The four different clas- 
sical phases are labeled by color: Blue - FM, green - NAF, 
orange - CAFa, red - CAFb. The thick lines correspond to 
the classical phase boundaries, the symbols indicate the loca- 
tions of the parameters used in Fig. and are labeled with the 
corresponding character, also used in Table [J The white dot 
represents the standard nearest-neighbor Heisenberg model 
(Jia = Jk, = Ji, J2 = 0), the black dots denote experimental 
points for BaFe2As2 in Ref . [y and CaFe2As2 in Refs. (|g| and 

m 



B(k) = —2 [Jib cos(k y b) + 2 J 2 cos(fc x a) cos(k y b)] , 

Figure Q] shows a contour plot of the total ground-state 
energy E gs = 1-L C \ + T-l Z p at zero field as a function of 
the frustration angle <f> and the anisotropy parameter 9. 
The energy unit is J c . Four magnetic phases appear (see 
caption) in a symmetric pattern in the <j>, 0-plane. We 
present the complete phase diagram of the Ji a ,6- J 2 model 
although only the sector < 9/n < 0.25, < <j>/ir < 0.5 
seems to be relevant for the Fe pnictide class according 
to Table fl] We notice the following characteristics of the 
phase diagram: 

i) The ground state energy and phase diagram are in- 
variant under the following symmetry transformations: 
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FIG. 2. Ground-state energy as a function of the frustration 
angle 4> f° r fixed anisotropy parameter 8 — n/8 from linear 
spin-wave theory (solid line) and exact diagonalization (solid 
dots; 24 sites). For comparison, the dashed curve shows the 
classical ground-state energy. Here and in subsequent plots of 
4> dependent quantities the vertical lines indicate the classical 
phase boundaries. 



System Ref. g SJ la SJ lb SJ 2 SJ C 4>/-k 8/tv 

1 CaFe 2 As 2 £ ■ 41 10 21 36 0.19 0.08 

C CaFe 2 As 2 J} 0.4 24-37 7-20 28-34 33-45 0.29 0.13 

D CaFe 2 As 2 W 0.22 49.9 -5.7 18.9 53.7 0.11 -0.04 

B BaFe 2 As 2 j3 0.28 17.5 17.5 35 39.1 0.35 0.25 

5 BaFe 2 As 2 j) 0.54 36 -7 18 31.6 0.19 -0.06 

6 CaFe 2 As 2 j> 0.75 27.4 -2.1 14.5 24.3 0.20 -0.02 

7 BaFe 2 As 2 5 0.84 36.1 -2.6 12.0 38.0 0.10 -0.02 

8 SrFe 2 As 2 £ 0.84 35.3 2.2 13.4 28.4 0.16 0.02 

TABLE I. Fe pnictide moment n — 25/iB and exchange inter- 
actions (in meV) from experiment (top) and theory (bottom) . 
Here J c is the average exchange energy scale and 8, <j> are 
anisotropy and frustration angles (Eq. l[5|). The first column 
holds the labels used in Fig. [T] (letters) and Fig. [15] (letters and 
numbers). 

tice with size N = 24 which is quite close to the spin 
wave results. 



Reflections at the lines 9 — ir/4 and — 37r/4 and inver- 
sion at the points (</>, 9) = (±7t/2, 3tt/4). Both opera- 
tions lead to (Jia,Jib) —> {Jib, J ia) with J 2 unchanged. 
This corresponds to an interchange of the the columnar 
CAFa/b phases while FM and NAF are mapped identi- 
cally. The classical ground state energy has even more 
symmetries. 

ii) In the isotropic case {9 — 7r/4) CAFa and CAFb 
are degenerate, and moving away from this symmetry 
line one of the two phases is selected. The stability re- 
gion of the columnar phases along the frustration axis 
{(/)) increases upon going away from the symmetry line 
9 = 7r/4 while that of the neighboring NAF or FM phase 
decreases. Therefore CAFa/b is stabilized by the pres- 
ence of a Ji a ,b anisotropy. 

iii) The exchange frustration is largest where three 
phases (e.g. the corner point 9/tt = 0.25, 4>/n ~ 0.15 and 
J2/J1 = 1/2) or two phases (the boundary lines) meet. 
Therefore the degree of frustration in a given compound 
in Table U also depends on the size of its anisotropy. 
While the above corner point is strongly frustrated and 
in fact has no long range order (see Sect. IIV C|) the point 
9 = 0, 0/tt as 0.15 {J lb = 0, close to (D) in Fig. [TJ is 
not strongly frustrated but stable CAFa despite having 
0/vr a 0.15 or J 2 /{Ji a /V2) = 1/2. 

Apart from the ferromagnet, which is an eigenstate to 
the full Hamiltonian, the spin-wave corrections stabilize 
the classical ground state, i.e., the zero-point energy in 
Eq. (j2"0"l) is negative semidefinite for all values of <j> and 
9. As an example, Fig. [5] displays the dependence of the 
ground-state energy on the frustration angle <fi obtained 
both from linear spin-wave theory (solid line) and from 
the classical model (dashed line). The plot was made 
with fixed anisotropy parameter 9 = 7r/8, corresponding 
to a spatial anisotropy Jn/Jia = v2 — 1 ~ 0.414. For 
a comparison we also present the numerical results from 
Lanczos calculations for finite clusters of the square lat- 



B. Spin- wave spectra 

We shall now discuss the excitation spectrum of the 
Hamiltonian (TK)1) for some typical points in the {</), 9) 
phase diagram, Fig. Q] In Table HI we have compiled 
an excerpt from the available literature on experimen- 
tal and theoretical values for the exchange parameters of 
j4Fe2As2 compounds, where A denotes an alkaline metal. 
From this list, we have chosen three parameter sets, indi- 
cated by the black dots in Fig. [IJ and the standard Neel 
antiferromagnet (white dot). 

Figure [3] shows plots of the k dependence of the spin 
wave excitations for different parameter sets {<fi,9). The 
parameter sets used for the plots are indicated by the 
symbols in Fig. Q] For simplicity, we scale k with the 
lattice constants and set k x a —> k x and k v b — > k y . All 
plots refer to the full crystallographic Brillouin zone. 

The top left spectrum (A) in Fig. [3] shows the well- 
known dispersion for the nearest-neighbor Heisenberg 
model for comparison. It has a Goldstone mode at the 
equivalent wave vectors Q — and Q = (±7r, ±7r). The 
low-energy dispersion u){k) = SE{k) is linear around 
these points with 

w(fc) = 25V Jia + Jib (28) 
x yj {J la - 2 J 2 ) {k x - Q x f + { J lb - 2J 2 ) {k y - Q y f. 

The top-right plot (B) shows the dispersion for 
{4>,9)/tt = (0.35,0.25), corresponding to an isotropic ex- 
change on the border between CAFa and CAFb phases. 
These parameters have been determined for BaFe2As2 
in Ref. 0. The Goldstone modes are at Q = (0, ±7r) 
and Q = (±7T, 0) and the equivalent points Q = and 
Q = (±7r,±7r), reflecting the twofold degeneracy of the 



FIG. 3. Spin wave spectra u(k)/J c . Clockwise from top left: (A) ((f), 9) /n — (0,0.25) - NAF, isotropic exchange with J2 = 0; 
(B) (<f>,0)/Tr = (0.35, 0.25) - CAFa/b, isotropic exchange; (C) (0, 6>)/tt = (0.29, 0.13) - inside CAFa; (D) (<£,0)/ir = (0.11, -0.04) 
- even more inside CAFa. 



CAFa and CAFb phases. The linear dispersion around 
the minima is 

tj a (k) =2S^/2J 2 + J la (29) 

x yj (2 J 2 + J la ) (k x - Q x f + (2 J 2 - J lb ) {k y - Q v f 

for the CAFa phase, and 

Wb(fc) = 2Sy/2J 2 + Ji b (30) 

x yj (2 J 2 - J la ) (k x - Q x f + (2J 2 + J lb ) [k y - Q y f 

for the CAFb phase. 

For (4>,0)/n = (0.29,0.13), assigned to CaFe 2 As 2 in 
Ref. [^, we show the spin- wave dispersion in the bottom 
right plot (C) of Fig. [3] With these parameters, the 
system is deep inside the CAFa phase. In contrast to 
the isotropic case, the dispersion around Q* = (0, ±7r) 
and (±7r, ±7r), while still being local minima (but with a 
quadratic k dependence), have a finite energy gap. We 
have E a (k) — remaining only at the wave vectors k = 
and (±7r, 0), characteristic for the CAFa phase. 

Finally, the bottom left plot (D) in Fig. [3] displays the 
dispersion for (cj),0)/n = (0.11,-0.04). This alternative 
parameter set was proposed in Ref. Hoi for CaFe 2 As 2 . The 



local minima at Q* — (0, ±7r) and (±7r,±7r) discussed 
in the previous paragraph have almost disappeared, the 
dispersion at the zone boundary k y = ±7r is flat. This 
property can be utilized to decide between the two ap- 
parently different parameter sets for the identical com- 
pound. In fact from this comparison it was concluded 10 
that the strongly anisotropic set (D) describes the dis- 
persion along {0,k y ) much better for large wavevectors 
ky/n > 0.5. Since the ordering is still of the CAFa type, 
the Goldstone mode at Q — (±7r, 0) remains for the 2D 
model. A realistic description of the dispersion requires, 
however, the inclusion of interplane exchange-' — which 
leads to a finite gap at these points^ 



C. Ordered moment 

The most appropriate quantity for judging the degree 
of frustration in the local moment model is the size of the 
ordered ground state moment m s (</>, 9) relative to its size 
for the unfrustrated (J 2 = 0) isotropic (J\ a = Jib = J\) 
NAF state. The latter is already reduced with respect to 
the classical value S to mj? ~ 0.6075 1 . The stronger the 
frustration the more m s should be reduced even below 
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the NAF value. In the isotropic case when J 2 < there 
is obviously no frustration and the NAF state is even 
stabilized. The essential question is: how large is the 
degree of frustration in the CAFb,b states relevant for 
the Fe pnictides? This question can be answered by cal- 
culating the moment reduction m s (<j), 0)/S in spin wave 
approximation. 

The ordered moment is the ground-state expectation 
value of the z component of the spin S' in the local co- 
ordinate system, 



1 

N 



S- 



1 

N 



E (4 a £ 



1 



k 



(31) 



Inserting the expression for unrequired to bring H into 
diagonal form (see Appendix |B] for a complete expres- 
sion), we get 



= S 



1 f ly A(k) - B(k) cos 2 8 C 



1 



E(h,k) 

' (32) 

Due to quantum fluctuations m s < S is smaller than 
in the classical case, except for the ferromagnet. Near 
the borders of the CAFb and CAFa phases to the NAF 
phase, the ordered moment vanishes, indicating the fail- 
ure of spin wave theory due to strong frustrations. Also 
between the FM and CAF phases the latter lead to a 
vanishing m s in a small region. 

Fig. [4] displays the behavior of m s = fi/fiB as a function 
of the frustration angle 4> for three different anisotropy 
parameters 9. The upper plot shows the isotropic case 
(see also Ref.[H), 9 = 7r/4, corresponding to J la — Ju = 
J\. Coming from the FM phase for 4> < —ir/2, m s is grad- 
ually suppressed to the well-known value m' « 0.3034 at 
= (J2 = 0), corresponding to the NAF state of the 
nearest-neighbor Heisenberg model. 

Increasing <f> towards the NAF-CAFa/b boundary fur- 
ther reduces m s , until at 4>/n ~ 0.1150 (where J2/J1 ~ 
0.3779), the ordered moment vanishes. The classical 
CAFa/b-NAF border is given by J2/J1 = 1/2, or 
(j>/Tr w 0.1476. Soon after entering the CAFa/b regime, 
at <p/ir s» 0.1508 or J2/J1 ~ 0.5129, m s becomes finite 
again and grows rapidly towards an almost constant value 
m s « m s (Ji = 0). 

At Ji = or = 7r/2, the lattice can be subdi- 
vided into two noninteracting sublattices with a nearest- 
neighbor interaction J = J2, therefore we must have 
m s (<fi = 7r/2) = m s ((f) = 0). This is indicated by the 
dashed horizontal line which illustrates that throughout 
most of the CAF region the moment reduction is almost 
the same as that of the unfrustrated NAF. In fact for 
(j) = 7r/2 the CAF moment is stabilized by quantum fluc- 
tuations which orient the moments of the two sublattices 
parallel. This is the so-called 'order by disorder' mecha- 




FIG. 4. Ordered moment m s = fi/fis for fixed anisotropy 
parameter 6 as a function of the frustration angle </)■ Top: 

8 — n j 4 (isotropic case, Ji a ~ Jit). The ordered moment 
vanishes for 0.1150 < (j>/n < 0.1508 (NAF-CAFa/b) and 
0.8491 < <j>/iv < 0.8524 (CAFa/b-FM crossover). Middle: 

9 — 7r/8. The gap between NAF and CAFa is smaller but 
still finite, the disordered region between CAFb and FM has 
disappeared. Bottom: 9 = 0, corresponding to Jk, = 0. The 
gap m s = appears for -0.00727 < 4>/-k < 0.00775. The 
dashed horizontal line in the three plots denotes the value 
m s ~ 0.3034 obtained for the standard Heisenberg model with 
Jia = Jib and J 2 = 0. 



The behavior of the ordered moment in the CAFa/b 
phase is independent of the sign of J\ and therefore sym- 
metric around (f> = ir/2. Reaching the border to the FM 
phase, m s sharply drops back to zero at <j>/-K f=s 0.8492 
or J2/J1 « -0.5129. The CAFa/b-FM border is given 
by J2/J1 = -1/2, or 0/tt rj 0.8524. At this border, m s 
immediately jumps to the saturation value m s = 1/2 in 
the FM phase. 

Now we turn to the anisotropic case, 9 ^ 7r/4. The 
lower two plots of Fig. @] show m s for 9 = 7r/8, corre- 
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FIG. 5. Ordered moment m s for cj>/ir = 0.852 as a function of 
the anisotropy parameter 0. The frustration angle is chosen 
such that in the isotropic case (9 = tt/4), the system is in the 
disordered regime at the CAFa/b-FM corner. 

sponding to J\bjJ\ a = v2 — 1 ~ 0.41, and the fully 
anisotropic case 9 — 0, meaning Jib = 0. For (jy/n < 1/2, 
the overall behavior of m s is similar to the isotropic case: 
After leaving the FM regime, m s is suppressed down to 
zero, and becomes finite again after the crossover to the 
CAFa phase. However, there are two quantitative differ- 
ences: Firstly, the region where m s = is smaller than 
for 9 = 7r/4. Secondly and most importantly inside the 
columnar AF phases, the ordered moment is restored to 
even larger values than for isotropic exchange. 

Exactly at <fi = 7r/2, we have J\ a = Jib = 0, and 
m s shows a dip with the universal value m s sa 0.3034 
of the nearest-neighbor Heisenberg model, indicated by 
the dashed lines in the plots. For any value of 9, at this 
point only Ji is finite, and the same argument as for 
the isotropic case applies; we must have m s ((/) = n/2) = 
m s (4> = 0) for arbitrary ratios Jibj J\ a - 

In contrast to the isotropic case, the symmetry around 
the point 4> = tt/2 is lost, and m s is restored in the CAFb 
phase towards the saturation value upon entering the FM 
phase. There is no region around the CAFb/FM border 
where m s is suppressed as in the isotropic case. 

A vanishing ordered moment implies that, at least 
within our approximation, the order parameter for the 
corresponding classical phase is destroyed by quantum 
fluctuations. Historically this was one of the first indica- 
tions of the appearance of an intermediate phase without 
magnetic order, and our findings suggest that the well- 
known disordered phase for the isotropic J1-J2 model for 
AF exchange couplings extends to the whole range of 
anisotropic interactions with arbitrary ratios J\b/ Jia- 

This is not the case for the CAFa/b-FM crossover, 
where numerically already at a deviation A9/9 < 0.01 
from the isotropic value 9 = it /A the ordered moment 
remains well-defined around the classical CAFa/b- 
FM transition point. Fig. [5] illustrates this behavior: 
The plot shows the ordered moment as a function of 
the anisotropy parameter for fixed frustration angle 
4>/tt = 0.8520 in the whole phase diagram. The two 
sharp dips at 9 = — Sit /A and 7r/4 correspond to the 



disordered regimes at the CAFa/b-NAF corner and the 
CAFa/b-FM corner in the phase diagram, respectively. 
At the CAFa/b-FM corner around 9 = tt/4, we have 
m s = only for a tiny range 0.2493 < 9/n < 0.2507. 
However, the precise relation to the extension of the dis- 
ordered phase around this corner for finite orthorhombic 
anisotropy remains unclear. 

In summary, if one regards the lower two panels in 
Fig0]representing the anisotropic case one observes a re- 
markable fact: The moment reduction by quantum fluc- 
tuation in the CAFa/b phases is less than in the unfrus- 
trated simple nearest-neighbor NAF phase (open circle), 
except for a very small region close to the strongly frus- 
trated CAFa/NAF boundary line. If one compiles the re- 
duced moments m s (<p, 9) for the proposed parameter sets 
of Fe-pnictide compounds (Table QJ as shown in Fig. [6] it 
is obvious that in most cases the moment reduction by 
quantum fluctuations for the proposed CAFa models is 
less than in the simple nearest-neighbor NAF 

This result is in part due to the stabilization of the 
moment due to the effect of the anisotropy as visible 
from Fig. [H which extends the stable range of 4> for the 
CAFa phase in Fig. [T] In fact the frustration angle for 
BaFe2As2 (D) is rather close to the strongly frustrated 
value </>/7r=0.15 of the isotropic (9/tt = 0.25) model; 
nevertheless it is at a considerable distance from the 
anisotropic CAFa/NAF instability line and hence shows 
only moderate moment reduction. Therefore from Fig. [6] 
we conclude that frustration/quantum fluctuation effects 
within a local moment picture may not be used to explain 
the surprisingly small ordered moment of the pnictides. 
We note, however that this conclusion does not invalidate 
the usefulness of the J\ a ,b~Ji local moment model for the 
interpretation of INS spin wave results. In classical (lin- 
ear) spin wave theory only the products SJi enter the 
spin wave velocity and dispersion SE(k) and therefore 
the shape of the dispersion does not depend on the size 
of the staggered moment as long as this approximation 
is reasonable. 



V. FIELD DEPENDENT PROPERTIES OF THE 
Jia, b -J2 MODEL 

It has been shown that high field magnetization (up 
to the saturation field h s ) is an excellent way to investi- 
gate the isotropic Ji- J 2 mode l 33 ' 34 because the degree of 
exchange frustration determines the characteristic non- 
linearity of the magnetization curve. This is especially 
true for the parameter region where the CAF phase is 
realized. One should remark however that in the un- 
doped Fe-pnictides the energy scale of J c ~ 5 • 10 2 K is 
too large to reach the high field region. However it might 
be feasible for some of the doped Fe-pnictides where the 
ordering temperature T m and hence J c are strongly re- 
duced, provided that doping in the spacing layers does 
not impede the usefulness of the J\ a .b-Ji model for the 
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- 0.4 - 



compound index 

FIG. 6. Ordered moments m s (<t>,8) normalized to the clas- 
sical constant m c \ = S for the compounds listed in Table U 
calculated with Eq. (|32l) . The dashed horizontal line indi- 
cates m B {4> — 0,8 — 7r/4) for the simple J2 = 0, Ji a = Jit 
nearest-neighbor Heisenberg model (Point A in Fig. [T]). 
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spin excitations in the FeAs layers. 



A. Magnetization at high fields 

The total magnetization of the system is the ground- 
state expectation value of the z component of the spin S 
in the global coordinate system, 



(33) 



Since this is just the projection of the ordered magnetic 
moments onto the field direction, we can also write 



= m s cos 0. 



(34) 



Here, O is not the classical canting angle C [this would 
describe the field dependence of the classical spin system, 
i.e., a straight line m(h) = S(h/h s )], but rather must 
include the first-order corrections from linear spinwave 
theory. We thus have to regard as independent variable 
again, return to the Hamiltonian given by Eq. (1191) before 
the replacement — > C and minimize its corresponding 
ground-state energy with respect to 0. 

Since l/S corrections are already included in the 
ground-state energy E gs (Q c ) = H c i + Hz P of the linear 
spin- wave Hamiltonian given by Eqs. (fT0|) . (p~9|) . pO)) . 
and (|23p . we can equivalcntly use the definition of the 
total magnetization per site as the negative field deriva- 
tive of E gs (Q c ), 



(35) 



with C given by Eq. ©. The result is 
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FIG. 7. Uniform magnetic moment m per site as a function 
of the applied magnetic field h/h s at three different frustra- 
tion angles in the CAF phases, (f>/iv = 0.16 (near NAF), 0.25 
(CAFa), and 0.65 (CAFb near FM). Between each pair of ad- 
jacent curves an offset Am = 0.1 is inserted. The solid lines 
denote the field dependence in the isotropic case, 9 = 7r/4, the 
dashed lines denote the maximally anisotropic case, 9 — 0. 



(see Appendix [Cl for details). 

Fig. [7] displays three sets of curves of m(h) for different 
frustration angles. The magnetic field is normalized to 
the respective saturation field. The solid curves show 
the field dependence of the induced moment m in the 
isotropic case — 7r/4, J\ a = Jif,, and the dashes curves 
show the same quantity in the maximally anisotropic case 
9 = or Jib = 0. 

Deep inside the ordered phases (well separated from 
phase boundaries), the first-order corrections to the to- 
tal moment are small, as indicated by the small bending 
of m(h) for </>/7r — 0.25 (middle curve) in Fig. [7] It is 
only near the crossover between adjacent different phases, 
where corrections become strong, and the field depen- 
dence of the magnetization differs significantly from the 
classical behavior, as shown in m(h) for ip/ir = 0.65 with 
isotropic exchange constants. 

Introducing an anisotropy generally reduces the quan- 
tum corrections leading to the nonlinear magnetization 
(dashed curves) . According to Table U the Fe pnictides 
are not close to the phase boundaries, therefore the effect 
of the anisotropy on the nonlinear magnetization curves 
is not very prominent as indeed suggested by Fig. [7J 



B. Magnetic susceptibility at low fields 

Another sensitive probe to the degree of frustration is 
the low field uniform susceptibility \ — dm/dh which 
may be used to obtain further insight as an alternative 
to the ordered moment. We obtain 
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(37) 




0.4 

% 0.3 



0.2 



0.1 



FM 




CAFa 


CAFbj 


9/ 71 = 









0.0 

<f>l7l 



FIG. 8. Uniform magnetic susceptibility \ — {^ m /9h) h as 
a function of the frustration angle cj> for (from top to bottom) 
9/n = 1/4, 9/n = 0.15, and 9 = 0. 



where the terms in brackets include the first-order cor- 
rections to the classical value, 



1 S 
Y r i = — — — = — = const. 
A 2A(0) h s 



(38) 



In the FM regions, x is undefined, since the system at 
T = is already fully polarized. In the three AF phases, 
X diverges near the FM phase, apart from those areas, 
where already the vanishing ordered moment discussed 
above indicates that the ordered phase is no longer stable. 

Fig. [5] shows the dependence of x 011 the frustration 
angle <f> for three different anisotropy parameters 6. In 
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FIG. 9. Temperature dependence of the uniform magnetic 
susceptibility for <j)/n = 0.35, and 9 — 0.25 as for (B) in 
Table U with J c = 58.5 meV, corresponding to a temperature 
variation between 170 K and 510 K. 



the top of Fig. [5J the isotropic case is shown. The sus- 
ceptibility diverges at the crossover from the NAF to the 
FM phase (J la = J lb = 0, 0/tt = -1/2). Around the 
border to the CAF phases (</> = tan _1 (l/2)), it vanishes 
and becomes undefined. Qualitatively the same happens 
near the border to the FM phase at <fr = it — tan _1 (l/2). 
Apart from shifting the classical phase boundaries, intro- 
ducing an orthorhombic anisotropy generally stabilizes 
the magnetic CAF ground state. Therefore the 'gaps' 
where strong frustration destroys the magnetic order are 
gradually closed, see the center plot of Fig. [5] 

At the CAFb/FM crossover, the behavior of x even 
changes to a divergence upon increasing the orthorhom- 
bic asymmetry by lowering 9. The bottom part of Fig. [5] 
shows x f° r Jib = 0, where the "gap" to the FM phase 
is closed. This is fully compatible with the rapid sta- 
bilization of the CAFb ordered moment as function of 
increasing anisotropy (see. Fig 0} at the boundary to 
FM. 

At 4> — tt/2) the <j> dependence of the magnetic sus- 
ceptibility has the same feature as the ordered moment 
discussed in Sec. lIVCl At0 = 7r/2, Ji = 0, and we must 
have x{4> = 0,0 = 7r/4) = x(0 = tt/2,0 = arbitrary), 
therefore a small dip appears in the orthorhombic case. 

In the CAFa sector relevant for the pnictides the sus- 
ceptibility has a plateau value except very close to the 
gap of instability. The value is almost equal to that for 
the unfrustrated simple n.n. NAF. This underlines again 
that quantum fluctuations due to frustration cannot ex- 
plain the anomalous magnetism of Fe pnictides. 

Furthermore the uniform susceptibility was found to 
have an unexpected temperature dependence. Within 
the (semi-)metallic itinerant model characterized by elec- 
tron and hole pockets one would naively expect a con- 
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stant Pauli susceptibility above the ordering temperature 
T m and a reduction below due to the gap opening. While 
the latter was found for numerous pnictide compounds, 
the susceptibility above T m is not constant but still in- 
creases roughly linearly with temperature 3 ^. Explana- 
tions for this observation were given within the nonin- 
teracting two band model 3 ^, an interacting Fermi liquid 
picture including nonanalytic correction terms 3 ^, within 
a model of coexisting itinerant and localized moments^. 

On the other hand INS results have shown that low 
energy spin excitations can be well described by a suit- 
ably parametrized J\ a j>-Ji local moment model accord- 
ing to Table |U Then one should expect something sim- 
ilar for the low (zero) frequency susceptibility, at least 
qualitatively. To check this conjecture we performed fi- 
nite temperature Lanczos calculations as, e.g., described 
in Ref. |35| for various finite J\ a .b-Ji clusters. The re- 
sult is shown in Fig. [9] for a parameter set (B) corre- 
sponding to BaFe2As2- For temperatures T < 0.25J C 
the calculation becomes unreliable due to finite size ef- 
fects. The temperature range corresponds to ~ 170 K at 
the lower and ~ 510 K at the upper boundary. The in- 
crease linearity in T observed for BaFe2As2 from 150 K to 
300 K is qualitatively reproduced although the absolute 
increase is too large. We mention that in the combined 
itinerant-localized model of Ref. the increased linear 
T-dependence was also attributed to the local moment 
contribution. At the very least our calculation shows 
clearly that in the present range of measurement one 
should not yet expect the high temperature Curie law 
x(T) ~ l/T for local moment systems. This should be 
expected only quite above the maximum temperature for 
X(T), which is at about 1.1 J c or 750 K in the case of 
Fig. EI 



VI. SUMMARY AND CONCLUSION 

The local moment model for Fe pnictides has been sur- 
prisingly useful to explain the low energy spin excitations 
obtained in INS phenomenologically, albeit with the as- 
sumption of possibly very anisotropic exchange. The lat- 
ter may have its microscopic origin in underlying orbital 
order as proposed i n 12 ' 26 but this is still unexplored. The 
very usefulness of the local moment picture may be a con- 
sequence of Hund's rule correlations in the multiorbital 
state of Fe pnictides-^. 

In this work we have investigated in detail the empir- 
ical localized moment J\ a .b~Ji model, in particular the 
effect of the in-plane anisotropy and the frustration ef- 
fect. It has been a recurrent topic to explain the com- 
paratively small ordered moments in Fe-pnictides as the 
effect of enhanced quantum fluctuations in the ground 
state due to large degeneracy caused by frustrated J\ a ,b 
and J 2 exchange bonds. 

We have investigated this question in detail us- 
ing spinwave approximation and in part the exact- 
diagonalization Lanczos method to calculate ground 



state energy, phase diagram and moment reduction by 
quantum fluctuation as function of anisotropy and frus- 
tration parameters. In addition we have studied high 
field magnetization and low field uniform susceptibility. 

We found that generally the anisotropy lifts the degen- 
eracy between CAFa/b phases and extends their stabil- 
ity range as a function of frustration. Furthermore the 
anisotropy reduces or closes the instability gap on the 
phase boundary to the NAF or FM phase respectively. 
Most importantly we have shown that in the CAFa sec- 
tor relevant for the pnictides according to Table U the 
moment reduction by quantum fluctuations is generally 
less than for the simple unfrustrated n.n. Neel antifer- 
romagnet. The same result can be obtained from the 
uniform low field susceptibility. 

Therefore we conclude that the anomalously low mo- 
ment in the pnictides is not explained by quantum fluctu- 
ations in effective localized moment models but needs a 
more microscopic viewpoint including the itinerant mul- 
tiorbital nature of the magnetic state. Such proposals 
have been made within recent ab-initio calculations us- 
ing the full orbital basis^r— . This does not invalidate, 
however, the exceptional usefulness of the simple J\ a ,b-J% 
model to describe the low energy spin excitations. 

We have also derived the spin wave excitations in an 
external field for the anisotropic model. It remains to be 
seen whether new information on the exchange models 
can be gained from INS experiments in finite fields. 

Finally, in a corollary we address a result of our anal- 
ysis not immediately relevant for pnictides because it is 
related to the magnetic instability at the CAFa/b-FM 
boundary {4>/it — 0.852). There are 2D local moment 
compounds 34 where the frustration angle is quite close 
to that boundary, contrary to the Fe pnictides. It has 
been shown for the isotropic model that the true ground 
state in this region is of the spin-nematic hidden order 
stated. Although spin wave theory is not adequate to 
fully address this question we have shown (Fig. [5]) that 
the columnar order at the boundary recovers immediately 
when turning on even a tiny anisotropy of n.n. exchange 
constants J\ a ,b- Since small anisotropics usually exist in 
such compounds we predict that the spin nematic state 
of the isotropic J\-J% model will be very hard to find in 
a real compound. 



Appendix A: Linear spin-wave analysis 

The formal procedures described here closely follow 
and generalize those discussed in Refs. |43j and |44|. where 
linear spin- wave theory has been applied to the triangular 
lattice. In the sections of this appendix, we don't make 
any assumptions about lattice geometry, dimensionality, 
and exchange constants except from the requirement that 
U(l) symmetry is conserved and still exists upon switch- 
ing on a magnetic field. The Hamiltonian is assumed to 
have the general form given by Eqs. ([1]) and @. 

Dropping the primes (working in the local coordinate 
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system), we use boson operators and a\ and write 

— S tt^Oj, 



S+ = V2S 1 - 



ST = V2Sa] 1 



2S 



1/2 



2Scii, 
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2Sa 



Keeping only terms up to bilinear order in the boson 
operators, we expand the scalar products in Eq. ((T|). The 
Hamiltonian up to bilinear order then reads 



s 

+ 2 



\{A a i + a * a ]) ( J ij cos(QRij) (cos 2 6 + l) + Jfj sin 2 9 
(0*0,3 + ajatj (j± cos(Q J R iJ ■) (cos 2 9 - l) + J* sin 2 9 
2 (aja, + a]aA fj+ cos(Q J R y ) sin 2 9 + J?, cos 2 9 
2i fajaj - ajatj sin((3-R l:) ) cos 9 
iV2S (a,i — a\ — dj + cij] Jij sia(QRij) sin 9 



+ V25 (a, + at + aj + at J (jf } - J± cos(Qi?y) ) cos 9 sin 9 
5 



a[ai cos 9 



^a, + a|J sin 9 



In the above sum, the contribution 



I [-iV25 (a 4 - a] - aj + at) sin(Qi? li ) sin 9 



I 

is antisymmetric in the variables i and j (site indices) 
and therefore vanishes upon summation. 

Inserting a Fourier representation of the spin wave op- 
erators aj into the above equation, replacing the sum over 
bonds (ij) with a sum over sites i plus their neighbors n 
gives after some rearrangement 



S COS 9^ Y ( J n ~ J n COs(QR r , 



a + aj 
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+\ cos e ( a l a k + H a % - l ) ■ 



to write 

n = 



Performing the sum over the neighbors n, together with 
Eqs. ©, ©, (HTJ), and (0 eventually leads to Eq. (ITol) . 
We define 

4 = (4°-s) 

JV5(5 + 1) (jx(Q) + A{0) cos 2 e) 

~Nh (2S + 1) cose + 1 £ atjfgag, 

( H x + H & H 2 \ 
\ H 2 Hi-H*)' 

A(k) - cos 2 6 (B(k) + 2A(ofj + | cos6, 

B(k) (i - cos 2 e c ) , 

C(k) cose, 



Hi 
#2 
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dropping the part linear in {at , a^ =0 }- Since the oper- 
ators aj: are bosons, their commutation relations can be 
written as ala z a? — 1 where a z is the symplectic unit 
matrix (which in our case is identical to the z Pauli spin 



matrix). Assume U% is the matrix which diagonalizes 
Hamiltonian. From the requirement that this transfor- 
mation respects the canonical commutation relations, it 
follows that Us must be symplectic, too, 



and we have 



Uta z U R = a z 



where D is the diagonal form of He. From this transfor- 
mation, we get the eigenvalues of 



—H 2 —Hi + H a 



by evaluating the characteristic polynomial x(H) = 
det (a z H? — E). This polynomial can always be written 
as 



X {E) = det (E - H a f - (Hi - H 2 ) (Hi + H 2 ) 



and the desired spin-wave dispersion can be immediately 
read off, 



E(h, k) = ^(Hi - H 2 ) (Hi + H 2 ) + H a 



A(k) - B(k) - 2,4(0) cos 2 6 + - cos 6 



A(k) + B(k) (1-2 cos 2 6) - 2A(0) cos 2 Q + - cos 6 



A(k) - cos 2 e (B(ic) + 2A(0)) + - cose 



1/2 



+ C(k) cose 



B(fc)(i-cos 2 e) +c(fc)cose, 



where we have taken the root with the positive sign only. Appendix B: Ordered moment 

Setting 6 — e c eventually leads to the dispersion (|23|) . 



Writing the Hamiltonian, Eq. (|X6|) , with the spin wave 
operators at, Eqs. ([21]) and ([22]) . and diagonalizing the 
LSW Hamiltonian in a more explicit way we get 



signB(fc) 



\ 



A(k) - cos 2 e ( B(k) + 2A(0)J + | cos e 
E(h, k) 
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\ 



1 ( A{k) - cos 2 9 \B(k) + 24(0) J + § cos 6 
E(h, k) 



Setting O = O c and inserting the coefficients i>g into 
Eq. (|3ip yields the expression (p£Ij) for the ordered mo- 
ment. 



at this level of approximation is a consequence of the 
zero-point fluctuations only. Integrating Eq. (|35[) be- 
tween h = and h — h s , we can write 



Appendix C: Uniform moment and susceptibility 

The uniform moment is given by Eq. p5|) . or equiva- 
lently 



rn 



1 dE s3 (Q c ) Scos0 c 
iV <9cos© c 9ft 



(CI) 



with the ground-state energy in linear spin- wave approx- 
imation 



E gs (Q c )=NS 2 (j x (Q)-A(0) 



cos 0, 



NSJ ± (Q) + ^J2 E ( h >V 



and E{h, k) given by Eq. (|23|) . We have 

M = - M( o>cose« + §£fM, 

a cos 6 C 2 ^ o cos Be 

fc 

d cos 6 r 1 



9ft 2SA(0) ' 

5 ' =-2cose c - 
acos0 r 



E(h,k) 



Inserting these expressions into Eq. (|C 1|) gives the desired 
result, Eq. (|36]l. 

The correction to the magnetization (the deviation 
from the classical behavior to c i = S'cosOc = S(h/h s )) 



AE gs = £ gs (ft = ft s ) - £ gs (ft = 0) = -iV / m(ft)dft 



(C2) 

As discussed elsewhere, the zero-point fluctuations re- 
duce £?g S (ft = 0), but do not reduce E gs (h s ). For the 
latter, the zero-point fluctuations vanish, again because 
the fully polarized state at h s is an eigenstate of H. Thus 
the classical energy difference is smaller than the energy 
difference including first-order corrections in Eq. (|C2j) . 



AEg - AE^ W 

= E*(h s ) - Ef a (0) ~ 
= £ g L s sw (0)-i^(0) 
= NSJ^(Q) 



2 ^ 



E^(h s )-E^(Q)) 



E{k) 



E 2p < 0, 



hence the corrections to the integrated magnetization 
must be negative. Assuming a monotonic behavior of 
m(ft), the same relation holds for the kernel in Eq. (|C2I) . 
and it follows that 

mLSw(ft) < m c i(ft), h < h s , 

consequently we can expect that the corrected magneti- 
zation curve lies below the classical one. 

For the susceptibility, we differentiate m once more 
with respect to the applied field, 

dm d cos O r 



<9cosO r dh 



and using 



J 



dm 



<9cosO c cos C 



B(k) (A(k) - B(k)) dE{h ji) 
d cos C 



S 



1 1 X ^ 

1 + 2SN ^ 



B 



A(0)E 2 (h,k) 
(k) (A(k) - B(k) 



A(0)E(h,k) 



B 



'SN 



\k) (A(k)-B(k) 
A(0)E 3 (h,k) 



We obtain Eq. J37]). 
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